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Abstract 

The image quality from Ground-Layer Adaptive Optics (GLAO) can be gradu- 
ally increased with decreased contiguous field of view. This trade-off is dependent 
on the vertical profile of the optical turbulence (C^ profiles). It is known that the ac- 
curacy of the vertical distribution measured by existing profiling techniques is 
currently quite uncertain for wide field performance predictions 4 to 20 arcmin- 
utes. With assumed uncertainties in measurements from Generalized-SCIDAR 
(GS), SODAR plus MASS we quantify the impact of this uncertainty on the trade- 
off between field of view and image quality for photometry of science targets at the 
resolution limit. We use a point spread function (PSF) model defined analytically 
in the spatial frequency domain to compute the relevant photometry figure of merit 
at infrared wavelengths. Statistics of this PSF analysis on a database of mea- 
surements are presented for Mt. Graham, Arizona and Dome C, Antarctica. This 
research is part of the activities of ForOT (3D Forecasting of Optical Turbulence 
above astronomical sites). 
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1. Introduction 



Characterization of the optical turbulence in the first few kilometres above the telescope is 



important for predicting the performance of Ground-Layer Adaptive Optics (GLAO) telescopes 
as a function of field of view diameter. Systems that have been proposed will correct visible 
or near- infrared science fields that are typically 4 arcminutes, and potentially up to 20 arcmin- 
utes in diameter and contiguous. There are several measurement techniques being advanced to 
provide statistics on the vertical distribution of the structure function coefficient C^{h) , and 
in this paper we explore the impact of a potential bias from generalized-SCIDAR and MASS 
measurements. The first of two sites we will investigate is a typical mid-latitude observatory 
site, Mount Graham (32.7 N, 109.87 W, 3200 meters), measured with generalized-SCIDAR. 
There are conifer trees at the summit with a height similar to the SCIDAR telescope's primary 
mirror, about 8 meters above the ground. The second is Do meC (75.1 S, 123.3 E, 3 260 meters). 



an Antarctic site with MASS and SO DAR measurements by lLawrence et al.l(|2004) and balloon 



measurements by lAgabi et al. 



(120061) 
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The GLAO PSF figure of merit that is of particular importance to wide field astronomy 
is radius of 50% encircled energy, computed at several points in the contiguous field of view 
and then averaged. It will be symbolized as EE50 here. EE50 is very closely related to 
the integration time to achieve some signal to noise ratio in background-limited point source 



photometry in the field (|Andersen et al.ll2006|) . a rather common science application for fields 



of view 4 to 20 arcminutes in diameter. Roughly, 



integration time oc EE50^. 



(1) 



We will compute EE50 starting with an analytically defined p hase Power Spectral 



(PSD ) for anisoplanatism and fitting error using established theory (|Jolissaint et al 



Density 



2006 



Tokovinin 



20041). Table [H lists the model parameters selected here. Computation from the analytic PSD is 
a fast method to discover the performance gradient of EE50{6), where 9 is the diameter of the 
field of view. 

The exact range of altitudes in the first few kilometres where bias has greatest impact de- 
pends on the basic GLAO system parameters, namely the diameter of the guide star asterism 
(also 9) whose signal is averaged and the effective pitch that is controlled by the ground conju- 
gated deformable mirror (A). The ratio hcz = A/6' defines the altitude below which any contri- 



bution to anisoplanatism is negligible. The term gray-zone (GZ) was coined (|Tokovinin 



2004) 



to identify the altitudes above hcz, where the contribution to anisoplanatism is not negligible 
(also known as partially corrected zone)ll| Fig{T]helps illustrate this in terms of performance in 
the focal plane. The plot shows the EE50 figure of merit as a function of the height of one layer 



of turbulence added to a typical, smooth profile. The layer co ptair 



'Looking at the approximate error transfer function in equation (8) ofjokovinin 
is the case. 



' of th e total turbulence 
(120041) one can see why this 
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strength of the smooth profile. FigH] shows that the largest performance gradient is at altitudes 
just above hcz- The gradient vanishes above ho = D/9, where D is the telescope diameter. 
In the following sections we will re-compute EE50{9) with estimated bias in the proportion of 
turbulence attributed to heights above or below hcz- 



2. Mount Graham and Dome C profile monitoring data 



The Mt. Graham G-SCIDAR measurements include 851 in High Vertical Resolution (HVR) 
mode and 991 1 in regular mode, both have been reduced to discretized turbulence strength Jj 
at height hi. Th ese were computed from the normalized co variance function of the irradiance 



fluctuations (see 

toC'M by 



Egner et al. 



2006L 



2007|) which are proportional to Jj, which are in turn related 



/ dhCl{h). 



(2) 



The intrinsic vertical resolution of SCIDAR is roughly given by 



0.78 



P 



X\h + hn 



(3) 



where p is the binary separation (35"), A is the wavelength of the scintillation signal (0.5/im), 
and hgs is the conjugation height of the generalized SCIDAR analysis plane (about —3500m). 
The regular mode resolution will represent free-atmosphere, above 1000 meters. The current 
HVR data set samples the scale height of the boundary-layer and provides data up to 1000 me- 
ters altitude. In a subsequent section we will describe how the ground-layer and free-atmosphere 
are reduced to form a composite statistical model. 

F or Dome C we will use 1701 MASS-i-SODAR profile monitoring measurements at Dome 



Cby 



Lawrence et al. 



(|2004h during the Antarctic winter of 2004. These data sample only two 



grid points between 30 and 1000 meters and do not sample any t urbulence below 3 meters 



2006h that give 



However, there exist balloon-borne micro-thermal measurements (|Agabi et al. 
us an estimate of the scale height and total strength of the ground-layer, and with this infor- 
mation we model the statistics of eight grid points fro m a height of zero to 2 00 meters. The 



turbulence measurements recorded by SODAR in the iLawrence et alj (|2004|) data we appro 



priate to a slab concentrated at 250 meters between the modelled ground layer and the lowest 
MASS measurement at 500 meters. 

For the Dome C altitudes from zero to 200 meters we define the following exponential 
model to 

ClQi) = Ae^-^l^^\ (4) 

Using Eqn.© it follows that 



J, = -Ah A (e^-'^ 



(5) 



We will choose the boundaries hi,^ in ^ Using a average, weighted by C^{h) 



hi 



-Ah A 



dh cl(h) h 



r 



(6) 



It has been observed with balloon measurements at Cerro Pachon (iTokovinin and Travouillon 



2006|) that the strength of ground-layer is governed primarily by the scale height. In our model 
we will make the scale height dictate the strength exclusively. A lognormal distribution of 
values of the scale height, /i^, while A = 740. x 10^^^ and is fixed, will give a lognormal 



distribution in seeing. 
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The Mt. Graham (MG) scenario has weaker overall seeing (median 0.74 arcseconds) than 
Dome C (DC, median 1.2 arcseconds). To illustrate the differences in the vertical distributions 
for these two sites we reduce the data to cumulative histograms of seeing in three slabs, shown in 
Figj2l The Dome C free atmosphere (right panel) and even upper ground-layer slab (middle)are 
quite calm. Though the left and middle panels of Fig|2] are not proof, the scale height of the 



MG turbulence is resolved by the HV-GS technique in another analysis (|Egner et al.ll2006l) to be 
between 100 to 250 meters. The DC scenario clearly has most turbulence concentrated between 
the telescope and 30 meters range (left panel Fig©. 



3. Reduction to composite profiles 

Since the measurements of the ground-layer and free-atmosphere at these sites is not si- 
multaneous, we must create composite profiles that would closely reproduce the PSF statis- 
tics as though we had computed them on a full set of Ji{hi) data, uninterrupted in h and 
sampled at the same time . To d o this we sort and combine the profiles of as described in 



Tokovinin and TravouillonI (|2006r) using the assumption of uncorrelated ground-layer and free- 
atmosphere seeing. We will briefly re-describe the process here in the context of our data. 

The Mt. Graham HVR will provide the ground-layer below 1000 meters and the regular 
SCIDAR measurements will provide the free-atmosphere above 1000 meters. Three groups 
of profiles in the ground-layer are identified using the sum of Jj. The first group are those 
profiles within 5% of the 25*'' percentile are combined in a simple average for J^. We call 
them the "good" case. The SO*'' and 75*'' percentile profiles area combined similarly and called 
"typical" and "bad". In each group the grid of hi is identical and hence remains unchanged 
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by the combining process. The same process is done for the free-atmosphere. The result is a 
reduction to three ground layer profiles and three free-atmosphere profiles, which together have 
nine permutations for composite profiles that can reproduce the PSF statistics as though we had 
computed them on all of the Ji{hi) data. 

For Dome C we sort and combine the MASS-i-SODAR profile monitoring measurements 
of the free-atmosphere above 200 meters in the same way we described for Mt. Graham. The 
ground-layer model does not need to be sorted; the choice of three scale heights Ha = [14, 9, 22] 
meters provide the median, first and last quartile of the integrated ground-layer. 

4. Resampling the Composite Profiles 

In all cases the shape of the composite profiles, whether averaged over time or defined by a 
function is smooth and well sampled by the grid of Ji{hi) defined so far. Hence, we are permit- 
ted to resample the the Ji{hi) grid for the GLAO PSF model, which is affected by the density 
of points in the gray-zone. We increase the number of grid points in the gray-zone until the PSF 
figure of merit has reached an asymptote. This is trivial for the ground-layer of Dome C, we can 
define the grid and then re-compute Ji{hi) with Eqn.® and Eqn.®. For the measurements 
of Mount Graham and the free-atmosphere of Dome C we divide several measured Ji{hi) grid 
into more numerous Jj{hj) using linear interpolation of the original discretized C^{h) data. 
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5. Predicted GLAO performance gradient 

The reduced composite C^{h) profiles for each site are input for the computation of field 
averaged radius of 50% encircled energy of PSFs at a wavelength of 1.25/im, outlined in ^ and 
symbolized EE50. The aim is to asses the impact on GLAO performance by potential biases in 
the measured vertical distribution of the turbulence strength. We have selected the performance 
EE50{9) metric to do this. Figgis a 3x3 multi-panel plot showing EE50{9) at Mt. Graham 
(red) and Dome C (blue). The thicker lines are the median values while the thinner ones are the 
first and last quartiles of the ordinate. 

Let us first consider the central column of plots to identify the fundamental differences 
between weak and strong free-atmosphere sites. In the upper one we see the Mt. Graham 
(red) EE50 gracefully increasing with 9, as the bottom of the gray-zone reaches into the 
boundary-layer turbulence 100 to 250 meters thick. For this top middle panel the actuator pitch 
of the DM was 0.5 meters and the Dome C scenario only very weakly affected by anisopla- 
natism, a consequence of an inadequate number of actuators for that site. In the central panel 
the pitch is 0.38 meters, which improves correction at Mt. Graham slightly in all conditions, 
and greatly improves Dome C for median or better conditions. The median and first quartile 
EE50{9) curves of Dome C and Mt. Graham have similar shape because the ground-layer pro- 
files at Mt. Graham have similar exponential shape. The bottom plot shows the potential gain 
for Dome C when the wavefront is controlled to a pitch of 0.1 meters. In the central column of 
plots, the important distinction between the two sites is that Dome C is always under- actuated 
with A = 0.5 and sometimes near the diffraction-limited EE50 with A = O.L Mt. Graham 
on the other hand has more high altitude turbulence and is always limited by anisoplanatism for 
these A. 
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Next, consider the columns of panels to the left and right of FigjS] showing uncertainties 



pertinent to fie 



Tokovinin et al 



dof view trade-offs in GLAO telescope design. As indicated in figure 4 in 



(|2005r) both MASS and SCIDAR measurements are believed to produce faithful 
total integrals of turbulence, however, the vertical distribution may be biased. The left column 
of plots in Fig|3]were computed from the Ji{hi) times 0.5 in the domain /i^^ < hi < 6km, 
the balance was conserved by putting turbulence in the lowest layer, below hcz- Likewise the 
the right column of plots is Ji{hi) times 1.5 in the domain hg^ < hi < 6km, with the balance 
conserved by removing turbulence from the lowest layer. The change from the central column 
of plots to the left or the right is the slope of the curves, germane to designing a field of view 
trade-off. The performance of a wide field survey can be expressed using the number of square 
arcminutes of sky that can be imaged to some limiting magnitude per unit time. For an the- 
oretical seeing-limited telescope this is of course proportional to 6*^. For a GLAO telescope 
with field of view 6 it will be roughly proportional to {6/ EE50{9))'^. EE50{6) in the middle 
row of Figl3](A = 0.38 meters) the slope of the median Mt. Graham EE50{9) in the domain 
10 < 6* < 20 arcminutes is about 45% less or more in the left or right panels. It is about =f15% 
for Dome C. In terms of integration time{9) oc EE50{9y in the domain 10 < ^ < 20 we find 
the slope is ±60% for Mt. Graham, ±30% for Dome C. In other words, at a mid-latitude site 
similar to Mt. Graham, the predicted survey coverage of the GLAO telescope could potentially 
be wrong by as much as 60%. 
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6. Summary 



The GLAO telescope scenario simulated here is a common design for wide field science 
demanding a contiguous field. The estimate of 50% uncertainty in the proportion of turbulence 
strength between the the corrected-zone and the gray-zone (in the first 6 km) is based on a 
comparison between MASS and SCIDAR and here we calculate an uncertainty of 60% in the 
slope function EE50{9). Dome C is truly a unique site, and more immune to the 50% uncer- 
tainty. However, if the true uncertainly is not simply multiplicative the uncertainty propagated 
to EE50{9) for Dome C might be similar to that of Mt. Graham. 

Acknowledgements. 



We would like to thank the authors of lLawrence et al.l (|2004|) for providing their SODAR-i-MASS 



data. This work has been funded by the Marie Curie Excellence Grant (ForOT)-MEXT-CT- 
2005-023878. 



References 

Agabi, A., E. Aristidi, M. Azouit, E. Fossat, F. Martin, T. Sadibekova, J. Vernin, and A. Ziad, 
2006: First Whole Atmosphere Nighttime Seeing Measurements at Dome C, Antarctica. 
PASP, 118, 344-348, doi:10.1086/498728, ^arXiv: astro-ph/05104T8 

Andersen, D. R., J. Stoesz, S. Morris, M. Lloyd-Hart, D. Crampton, T. Butterley, B. EUerbroek, 
L. Jolissaint, N. M. Milton, R. Myers, K. Szeto, A. Tokovinin, J.-R Veran, and R. Wilson, 



10 



2006: Performance Modeling of a Wide-Field Ground-Layer Adaptive Optics System. PASP, 
118, 1574-1590, doi:10.1086/509266, [a7)riv:astro-ph/0 6100 97 

Egner, S. E., E. Masciadri, and D. McKenna, 2007: Generalized SCIDAR measurements at Mt. 
Graham. PASP accepted. 

Egner, S. E., E. Masciadri, D. McKenna, and T. M. Herbst, 2006: Beyond conventional G- 
SCIDAR: the ground-layer in high vertical resolution. Advances in Adaptive Optics II. Edited 
by Ellerbroek, Brent L.; Bonaccini Calia, Domenico. Proceedings of the SPIE, Volume 6272, 
pp. 627256 (2006)., Presented at the Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference, Vol. 6272, doi: 10.1 117/12.671380. 

Jolissaint, L., J. -P. Veran, and R. Conan, 2006: Analytical modeling of adaptive optics: foun- 
dations of the phase spatial power spectrum approach. Optical Society of America Journal A, 
23, 382-394. 

Lawrence, J. S., M. C. B. Ashley, A. Tokovinin, and T. Travouillon, 2004: Exceptional as- 
tronomical seeing conditions above Dome C in Antarctica. , 431, 278-281, doi: 10. 1038/ 
nature02929. 

Tokovinin, A., 2004: Seeing Improvement with Ground-Layer Adaptive Optics. PASP, 116, 
941-951, doi: 10. 1086/424805. 

Tokovinin, A. and T. Travouillon, 2006: Model of optical turbulence profile at Cerro Pachon. 
MNRAS, 365, 1235-1242, doi:10.1111/j.l365-2966.2005.09813.x. 

Tokovinin, A., J. Vernin, A. Ziad, and M. Chun, 2005: Optical Turbulence Profiles at Mauna 
Kea Measured by MASS and SCIDAR. PASP, 111, 395-400, doi: 10. 1086/428930. 



11 



of Figures 



The gray-zone begins above hcz 

Comparison of the Dome C (DC) and Mount Graham (MG) turbulence profile 

data used here 

The field averaged radius of 50% encircled energy on PSFs at 1.25/im, plotted 
as a function of the GLAO field of view 
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Fig. 2. Comparison of the Dome C (DC) and Mount Graham (MG) turbulence profile data 
used here. 
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Fig. 3. The field averaged radius of 50% encircled energy on PSFs at 1.25/im, plotted as a 
function of the GLAO field of view. 
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List of Tables 

1 The parameters and implicit assumptions of the GLAO PSF model 
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Table 1. The parameters and implicit assumptions of the GLAO PSF model, 

phase PSD von Karman, Lq = 30 meters 

telescope diameter D = 8 meters 

Beacons 4 point sources at range if = 90 km at zenith 

Beacons evenly distributed on a circle of diameter 9 in the field 

image wavelength A = 1.25//m 

image locations sampling a square field of view with vertices that intersect the circle 

Deformable Mirror cartesian grid of actuators with pitch, A 

Deformable Mirror each actuator has a sinc-like influence function 

Deformable Mirror conjugated to height = 
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